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A  Class  of  Energy  Levels  for  the  Heisenberg  Linear  Chain.  I 


Robert  B.  Griffiths* 


Abstract 

The  highest  and  lowest  energies  as  a  function  of  the 
total  spin  are  computed  for  the  claB3  of  "unbound"  states 
in  the  Be  the  formalism  for  the  linear  chain  of  spin  ^  atoms 
with  a  Heisenberg  exchange  interaction  between  nearest 
neighbors.  In  addition,  the  magnetization  as  a  function 
of  magnetic  field  is  calculated  in  the  limit  of  zero  tempera 
ture  for  an  infinite  antiferromagnetic  chain. 
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I.  Introduction 

The  properties  of  antiferromagnetic  insulators  are  often  discussed 
on  the  basis  of  the  Heisenberg  model  of  exchange  between  neighboring 
atoms.  Calculations  for  two  and  three  dimensional,  lattices  invariably 
proceed  by  means  of  approximations  whose  validity  is  difficult  to  Judge. 
Hence  there  is  still  some  interest  in  examining  the  one  dimensional  case, 
for  which  a  certain  amount  of  progress  has  been  made  toward  an  exact 
solution. 

Bethe^  showed  that  the  eigenvalue  problem  for  a  chain  of  N  spin  \ 
atoms  with  Hamiltonian 

*  «  17  i 

1*1 

•  5 ,  (l) 

where  is  the  spin  operator  for  the  l'th  atom  and  J  the  "exchange 
Integral,"  could  be  reduced  to  that  of  finding  the  solutions  to  a  set  of 
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coupled  transcendental  algebraic  equations.  Using  this  procedure,  Hulthdn 
calculated  the  exact  ground  state  energy  for  an  infinite  antiferromagnetic 
(J  >  0)  chain;  and  des  Cloizeaux  and  Pearson ^  have  recently  calculated 
the  energies  of  the  lowest -lying  excitations  or  "spin  waves." 

A  certain  subset  of  the  eigenstates  in  the  Be the  formalism  are 
"unbound";  this  nomenclature  is  made  more  precise  in  Section  II  below. 

The  unbound  states  are  of  particular  importance  in  discussing  the  proper¬ 
ties  of  the  antiferromagnetlc  chain,  since  they  include  the  antiferro¬ 
magnetic  ground  state  and  the  des  Cloizeaux  and  Pearson  "spin  waves." 


The  present  report  is  concerned  entirely  with  the  unbound  states.  In 
particular  ve  compute  In  Sections  III  B,  C,  and  D  the  minimum  and  maximum 
energies  for  unbound  states  having  a  given  total  spin  S,  as  a  function  of 
S.  The  minimum  energies  are  of  particular  interest  since  there  is  good 
reason  to  believe  that  they  also  represent  the  minimum  energies  of  all 
states  In  the  chain  (for  a  given  S).  This  permits  a  calculation  (in 
Section  IV)  of  the  magnetization  of  the  antiferromagnetic  chain  as  a  func¬ 
tion  of  magnetic  field  at  zero  temperature. 

The  Bethe  formalism  for  treating  the  eigenvalue  problem  is  summarized 
and  discussed  in  Section  II.  A  small  contribution  towards  the  formidable 
task  of  making  the  formalism  fully  rigorous  is  contained  in  an  existence 
proof  in  Appendix  A.  The  derivation  of  the  Hulthdn  Integral  equation, 
upon  which  the  work  in  Section  III  B,  C,  and  D  depends,  1b  summarized  in 
Section  III  A. 

The  results  of  our  computations  do  not  agree  with  a  previous  calcula¬ 
te 

tlon  by  Ledlnegg  and  Urban.  The  reason  for  the  discrepancy,  which  we 
believe  to  be  an  error  in  their  computations,  is  discussed  in  Appendix  B. 

By  considering  the  detailed  energy  level  structure  of  the  class  of 
"unbound"  states,  one  can  also  say  something  about  the  free  energy  of  the 
antiferromagnetic  chain  at  low  temperatures.  This  discussion  will  be 
deferred  to  another  report. 
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II.  The  Bethe  Equations 

Consider  a  linear  chain  of  N  spin  £  atoms,  closed  upon  itself,  with 
Hamiltonian  (l).  Assume  that  J  is  positive,  that  is,  the  interaction  1b 
antiferromagnetic.  The  largest  eigenvalue  of  ,  corresponding  to  the 
ferromagnetic  ground  state,  is 

5p  *  J  a/ct  ,  (2) 

while  the  lowest  eigenvalue,  the  energy  of  the  antiferromagnetic  ground 
state,  is2 

*  NJli-lhjl)  (3) 
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in  the  limit  of  very  large  N. 

Choose  as  basis  vectors  states  for  which  all  spins  In  the  chain  are 
either  "up"  or  "down"  with  respect  to  the  z  axis.  Let  f 

nln2* • ,nr 

be  a  state  for  which  the  spins  n^,  n^,  .  .  . ,  n^  are  down  and  all  other 
spins  are  up.  Any  state  with  z  component  of  spin  equal  to  N/2  -  r 
may  be  written  as  a  linear  combination  of  such  states: 


^  “  2-  a  (  •  ••  /  *r)  *».••• 


•V 


(*> 


where  the  summation  extends  over  all  sets  of  r  (distinct)  indices  n.. 

w 

The  eigenfunctions  discussed  by  Bethe1  are  of  the  form: 


K 


*»,  —‘-jflr)1 


r. 

I  e*p 
p-« 


'•I 


Y'i 


J<  4 


Kj ,  n ) 


(5) 
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where  the  summation  extends  over  all  r.'  permutations  of  the  integers 
1,  2,  .  .  . ,  r  among  themselves;  PJ  is  the  image  of  ,}  under  the  P'th 
permutation.  The  "vavevectors"  kj  satisfy  the  equations : 


Nfej-l^Aj+Z  }<f>u  j-  ,  r  («) 

where  the  ^j's  are  integers  between  0  and  N-l,  and  the  ^'s  are  defined  by: 


i  hi  m  i  i  i  hi) 

-*  *  -  « 


Let  E  be  the  energy  of  the  eigenstate,  and  define  the  "normalized 


energy"  e  by: 


g.  (e,  -£]/dThf), 


a  quantity  which  is  obviously  non -negative.  Similarly  define  by: 


€,-/»,  i  -  e  -  (£-£*,)/ (ijw). 


The  energy  of  the  state  (5)  is  given  by 


C./v-'Z.  ('-<•«*;) 


With  each  eigenstate  of  the  form  (5)  there  is  associated  a  set  of 
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Integers  Aj  through  (6).  The  order  of  the  A j  Is  obviously  unimportant. 
Furthermore,  a  state  with  sane  of  the  Xj  equal  to  zero  has  the  identical 
energy  and  total  spin  (though  not  the  same  z  component  of  spin)  as  the 
corresponding  state  for  which  all  the  zero  A 's  have  been  eliminated. ^ 

Thus  without  loss  of  generality  we  can  restrict  the  discussion  to  sets 
of  integers  such  that 

0  <A,  £  At  <  ..  .  £  A  r  <  N  .  (n) 

Since  the  chain  is  closed  upon  itself,  each  eigenstate  may  be  charac¬ 
terized  by  a  wave  number  q  given  by: 

<1  =  AT1  2  Aj  H)  =*z  hi  (mJ  in) ,  ^ 

j  i 

(Note  that  ^  =  -^.) 

Although  the  sum  of  all  the  k  is  real,  individual  k  values  may 

d 

be  real  or  complex.  This  fact,  among  others,  makes  the  discussion  of  the 
existence  of  roots  of  Eq.  (6)  rather  difficult.  The  following  questions 
have  not  been  discussed  (to  the  best  of  our  knowledge)  in  any  adequate 
sense:  For  what  choices  of  sets  of  integers  A,-  does  (6)  possess  a  solu¬ 
tion?  If  there  is  a  solution,  is  it  unique?  Given  a  solution  to  (6), 
what  are  the  conditions  such  that  the  wave  function  (5)  does  not  vanish? 
Are  the  wave  functions  corresponding  to  different  sets  of  Xj  linearly 
independent?  Can  all,  or  at  least  almost  all,  the  eigenfunctions  of  the 
chain  be  written  in  the  form  (5)?^ 
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It  Is  convenient  to  divide  the  eigenstates  of  the  form  (5)  into 
two  categories,  the  "bound"  states  and  the  "unbound"  states.  The  latter 
category  contains  the  states  for  which  all  the  k^  are  real,  the  former 
the  states  for  which  at  least  some  of  the  kj  are  complex.  There  is  a 
corresponding  division  of  the  sets  of  integers  { Aj  J  .  The  reason  for 
the  terminology  "bound"  and  "unbound"  is  to  be  found  in  Be the 's  discussion^ 
of  the  case  r  =  2.  He  shows  that  complex  k  values  correspond  to  a  state 
in  which  |a(n^,n2)j  decreases  exponentially  as  In^-n^l  increases;  i.e., 
two  spin  waves  interact  with  each  other  to  form  a  bound  state. 

Q 

In  particular  Bethe's  calculations  indicate  that  the  class  C°  con¬ 
sisting  of  those  sets  of  integers  for  which 

Aj4.,  z  \j  +- 1  (13) 

is  satisfied  are  unbound.  We  shall  denote  by  B  the  class  of  all  states 
for  which  (13)  is  not  satisfied.  Class  B  also  contains  some  unbound 
states,  although  Bethe's  computations  suggest  that  the  relative  number 
of  such  states  is  small. 

In  our  calculations  below  we  shall  be  concerned  entirely  with  states 
of  class  C.  We  shall  assume  that  each  set  of  integers  satisfying  (11)  and 
(13)  gives  rise  to  a  single  state  (or,  rather,  a  multlplet  of  degeneracy 
2S+1,  where  S  is  the  total  spin  of  the  state)  for  which  an  the  k  values 
axe  real,  and  that  states  corresponding  to  different  sets  of  Integers  in 
class  C  are  linearly  independent.  A  certain  (but  far  from  complete ) 
justification  for  these  assumptions  is  to  be  found  in  the  calculations  of 
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des  Cloizeaux  and  Pearson^  and  the  author^  on  finite  chains;  and  in  a 
proof  of  the  existence  of  solutions  to  Eq.  (6)  which  will  be  found  in 
Appendix  A. 

The  importance  of  the  states  in  class  C  for  the  problem  of  the  anti- 

ferromagnetic  chain  comes  from  the  following  "theorem":  For  a  given 

total  spin  S,  the  state  with  the  lowest  energy  (for  antiferromagnetic 

coupling)  belongs  to  class  C.  The  "theorem"  has  not  been  proved,  although 

it  is  supported  by  the  aforementioned  calculations  in  references  1  and  3, 

and  also  by  the  results  of  this  paper  (see  Section  III  D  and  Fig.  1). 

Some  insight  into  the  structure  of  Eq.  (6),  at  least  for  real  values 

of  the  k  ,  is  provided  by  regarding  it  as  a  non-linear  transformation 
J 

of  the  vector  (k^  ...»  kr)  into  another  vector  with  components: 


*/“  ZvN-'X;  +  VZ  4Ckj 


(14) 


where  for  S.,  we  have  written  A(k.,k  ).  A  solution  to  (6)  is  a  fixed  point 

T  J  JL 

of  the  transformation  (l4).  Suppose  that  we  attempt  to  solve  (14)  by- 
iteration.  A  sensible  starting  value  for  kj  is: 


kj,  *  Zrt/"1  Xj 


(15) 


which  lies  within  the  interval  (0,  2tt).  If  k^  and  k^  are  both  in  the 

y 

Hence  the  second  term  on  the  right  hand  side  of  (14)  represents  an  "attrac¬ 
tive  force"  between  pairs  of  vavevectors.  Thus  if,  for  example,  the  initial 


k^)  is  positive. 


interval  (0,  2rr)  and  k^  is  greater  than  k^,  then  q>(k 
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values  kjo  are  distributed  uniformly  over  the  interval  (0,  2 tt),  suc¬ 
cessive  iterations  of  (l4)  will  lead  to  an  increase  in  the  density  of 
the  kj ' s  near  tt  and  a  decrease  at  the  end  points  of  the  interval,  0  and 

2n. 


If  the  starting  values  for  k^  and  k^+1  are  too  close  together, 

successive  iterates  of  (14)  may  eventually  lead  to  a  "collision"  in  which 

kj+1  =  kj.  But  a  solution  of  (6)  for  which  two  of  the  k^'s  coincide  is 

trivial  in  the  sense  that  the  corresponding  wave function  (5)  vanishes 

Identically.  A  possible  remedy  for  this  situation  is  to  make  the  two 

kj's  complex,  with  imaginary  parts  of  opposite  sign,  resulting  in  a 

"bound  state."  Condition  (13)  for  states  in  class  C  insures  that  the 

initial  values  of  adjacent  k  are  far  enough  apart  so  that  no  "colli- 

J 

sions"  will  occur.  These  considerations  are  further  developed  in 
Appendix  A. 

As  a  first  step  in  examining  the  energy  level  structure  of  class  C 
we  shall  investigate  the  highest  and  lowest  energies  for  a  given  value 
of  S,  the  total  spin.  If  the  value  of  r  (and  thus  S)  is  fixed,  Eq.  (10) 
indicates  that  e  will  be  a  maximum  when  the  k^  are  near  the  center  of 
the  interval  (0,  2rr)  and  a  minimum  when  they  are  distributed  at  the  two 
ends  of  the  interval.  On  the  basis  of  (15)  and  the  discussion  following 
that  equation,  it  is  reasonable  to  assume  (although  we  have  no  proof) 
that  the  former  condition  is  achieved  when  the  Ay  are  near  the  center  of 
the  interval  (0,H)  and  the  latter  when  they  are  distributed  at  the  two 
ends  of  the  interval.  Further  recall  that  e  measures  the  energy 
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from  the  ferromagnetic  ground  state,  where as  we  are  Interested  In  neuimum 
and  minimum  energies  for  antiferromagnetic  coupling.  The  foregoing  argu¬ 
ment  leads  to  the  following  choices  for  a  given  value  of  S  =  ^N-r : 

Minimum  energy: 

A, «  iv-r a,,  w-,  Ar--iw*r-|.  (a6) 

In  particular,  the  ant i f e r romagne tic  ground  state  (S=0)  corresponds  to: 

A,  =  I ;  K  -  3 ,  K'C,  •  • .  j  *  m-\  .  (17) 

Maximum  energy: 

a  •  •  J  *r/t  *  r-lJ 

J  Xr/^+t*  U-ir-3)  t  J  Ar  m  ht-l  .  (l8) 

We  have  tacitly  assumed  In  writing  (l6),  (17),  and  (lfi)  that  N  and 
r  are  even;  minor  modifications  sure  required  if  one  or  the  other  Is  odd. 
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III.  The  Hulthtfn  Integral  Equation 

A.  The  Antiferromagnetic  Ground  State 

Consider  the  anti ferromagnetic  ground  state  which  corresponds  to 

the  set  of  Integers  (17).  For  large  N  it  Is  reasonable  to  replace  Eq.  (6) 

2 

by  an  integral  eolation  In  the  following  manner:  Let  x  be  a  continuous 
variable  on  the  interval  (0,1) .  Replace  A;  by  Nx,  by  k(x),  and 
Eq.  (6)  by 

N  kb)  -  lx  Nfi 

where 

cof  4  * 

*  Tt 


2.  4ihi)J y 


(19) 


4  [*« 4-  4  kl*)  -  ~ 

*4  <  k  . 


(20) 


The  normalized  energy  e  (see  Eq.  (10))  is  given  by: 


*  i  f  [l-  COS  (21) 

In  place  of  the  non-linear  Eq.  (19),  Hulthdn  obtained  a  linear  equa¬ 
tion  by  first  differentiating  both  sides  of  (19)  with  respect  to  x  and 
then  introducing  dx/dk  as  the  new  unknown  function  and  k  as  the  new 
independent  variable,  resulting  in  the  equation : 


4  +(o»t‘\,k  ,k') 


flnx{k 


(22) 
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Still  another  change  of  varlaoles: 

X  »  {.  h. 

-**U* 


permits  (22)  to  be  rewritten  in  the  form: 


Mi)  *  3*1*)  -  [ 

*« 

where 

/+ i*r' 


(23) 


(24) 


(25) 

(26) 


and,  furthermore, 

e.  *  Jf  ^ 


(27) 


HultbiSn  solved  (24)  by  successive  iteration.  For  our  purposes  it  is 


convenient  to  write  the  solution  in  the  form 

j.(r)  -J  *(■»-*»)  i.t'iuv 


(26) 


where  the  resolvent  kernel  R(x)  may  be  written  as  an  infinite  sum: 
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R(*)=  -L  t  »* 

>X  *•/  (?n)1  +  x'*  ^9) 

or  as  an  Integral: 

R(»)-f  (so 

i+  (i+x)* 

The  following  Is  an  important  relationship  between  the  kernel  K  and.  its 
resolvent  R : 

/  H(%-v)R(i-x)J'i*k(j-x)-Rlf-T;)  (31) 

The  solution  to  (24)  is: 

-  -J,  Sec.li  ■£  TT”^  (32) 

and  the  corresponding  energy  is: 

£0  *  !•]  1  .  (33) 


Ilf 


B.  States  of  Minimum  Energy 

Hultbrfn's  procedure  outlined  above  (with  some  changes  in  notation) 
may  be  applied  with  only  minor  modifications  to  the  more  general  problem 
of  finding  the  energies  associated  with  states  corresponding  to  the  sets 
of  X  values  in  (l6)  and  (l8).  Consider  first  the  states  of  minimum 
energy,  (l6).  Equations  (19)  and  (21)  must  be  replaced  by: 

f  t  +* 

k(*)  *  2  a  X  +-  {  J  (3|, 


u. " 


-  cos  k  C* )  7  ^  x 


with  a  determined  by  the  requirement  that 


-  S  =  r  =  Na 


In  place  of  r  and  S  it  is  convenient  to  introduce  the  abbreviations: 


r/N 

ft  -  S/b!  -  -p. 


Equation  (34)  may  be  transformed  by  the  same  procedure  used  with  (19). 


Define 


7  *  co+  ^4 
/frj  *  -dt/d'i  . 
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The  resulting  linear  integeral  equation  Is  : 


ft's)* 

-d 

with  limits  of  Integration  determined  by  the  requirement: 

riiA  ,<  r< 

r-<^ *  1J  *  *-{J  (i*U  TV?  f  HD# 

%-•  -<  J-d 

The  normalized  energy  s  is  equal  to: 

e- J*  f(f)  O+vr'-'t 


(39) 


(40) 


(M) 


The  quantity  of  interest  is,  of  course,  e  regarded  as  a  function  of  p  . 
The  most  convenient  procedure  to  follow  in  deducing  this  function, 
however,  is  to  choose  a  value  of  a  ,  solve  (39),  and  then  obtain  p 
and  e  by  means  of  (40)  and  (4l). 

In  Eqs.  (39)“(4l),  only  the  valueB  of  f (§)  for  |§|  <;  or  are  employed. 
However,  f (§)  may  be  defined  for  values  of  |g|  >  a  by  means  of  (39)  and 
thus  extended  to  a  function  on  the  entire  real  axis.  This  extended  func¬ 
tion  also  satisfies  a  different  integral  equation.  Multiply  the  rlgit 
and  left  hand  side  of  (39)  by  R(£  -  g)  and  integrate  with  respect  to  g 
from  -  oo  to  +  oo.  By  interchanging  the  order  of  Integration  on  the  right 
hand  side  and  utilizing  the  relations  (31),  (28},  and  (39);  one  obtains  the 
equation : 

j  f  f .  )  fh)Jy 

-Q*  J  tC  / 


(39a) 
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Similarly,  In  place  of  (40) 

1  * 

~  6  -  /T 

Equation  (40a)  Is  obtained  by  Integrating  both  sides  of  (39)  from  •  »  to 
+  <»,  noting  that 

!>>*-'  ^  ■  /  (tt> 

To  obtain  (4la),  multiply  the  right  and  left  hand  sides  of  (39a)  by 

(ihT1--  {rg.C-iJ 

and  Integrate  from  to  +  ».  By  Interchanging  the  order  of  Integration 

on  the  right  hand  side  and  utilizing  (27),  (28),  and  (4l),  one  arrives  at 
the  result  (4la). 

The  equations  (39a)-(4la)  are  particularly  useful  in  the  problem  of 
obtaining  the  asymptotic  behavior  of  and  when  a  is  large;  i.e. , 
when  «1  and  p1  are  small.  This  asymptotic  behavior  is  of  particular 
interest  since  it  determines  the  magnetic  susceptibility  in  small  magnetic 
fields  when  the  tenperatuze  is  zero,  as  shown  below  in  Section  IV.  An 
Investigation  of  this  problem  by  Ledinegg  and  Urban*  utilising  Eqp.  (39)- 
(4l)  appears  to  be  in  error  due  to  sense  unjustified  approximations.  We 
(  shall  discuss  their  calculation  further  in  Appendix  B. 


and  (41)  one  has: 


(40a) 

(43a) 


<  » --CM ri— 


IT 


Since  f(§)“f(-§)>  (39*)  may  be  rewritten  as: 

f(t)  ‘  f.(%)  +  O m 

Furthermore,  since 


<f.l i)  < e 


(<*) 


we  are  clearly  Justified,  in  replacing  fQ(§)  by  exp(-^rr§)  in  (1+3)  and 
(4la)  in  order  to  obtain  an  asymptotic  solution  for  large  a-  Mote  that 
only  value 8  of  5  and  T)  greater  than  a  come  into  consideration  in 
solving  the  integral  equation  and  evaluating  e-^  and  p  .  For  large 
values  of  z. 


R(*)»  R(-*)  /v  (ztk 


(*5) 


therefore,  provided  a  is  very  large,  it  should  he  possible  to  neglect 
R(5  +  T| )  relative  to  R(£  -  T|)  in  (1+3). 

Next  define: 


T>IV  -  e*”  f(r+«) 


<«) 


and  rewrite  (1+3),  using  the  approximations  Introduced  in  the  preceding 
paragraph,  in  the  form: 


pfa)  •  e'***  +J'K(*-v)i0 (y)df 


IS 


The  quantities  g^  end  are  given  by: 

ft  *  £  %1T‘l  (i  jo  V )*  «C  %T*  (48) 

e~T<(Til  C'*f  (*)<*/)  =  ie'^  •  (49) 


Numerical  solutions  of  (39)  confirm  the  exponential  variation  of 
and  predicted  by  (48)  and  (49),  with  coefficients  which  are 
slovly  varying  functions  of  a.  A  numerical  solution  of  (47),  though 
not  to  a  very  high  accuracy,  yields  the  following  values  for  a  and  b: 


a  =  .484  ±  .0005 
b  -  1.156  ±  .001 


(50) 


On  the  basis  of  a  model  ve  plan  to  discuss  in  a  subsequent  report,  we 
2  2 

expect  the  ratio  b/a  to  be  equal  to  ,  and  this  is  in  agreement  with 
(50)  within  the  precision  of  our  calculations.  Hence  we  write: 

e,  ~  \  **  V*2,  (51) 

for  very  small  p^. 

It  is  also  useful  to  consider  the  states  of  minimum  energy  for  the 
case  where  or  is  small,  and  hence  p  and  «  are  both  small.  Equations 
(39)-(4l)  yield  immediately  the  asymptotic  limits: 
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e  ~  **/•* 

€  ~  le  (52) 

provided  p  or  a  is  very  snail. 
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C.  States  of  Maximum  Energy 

For  a  given  value  of  total  spin  S,.  we  expect  the  states  with  the 
maximum  energy  in  class  C  to  be  associated,  with  sets  of  integers  of  the 
form  (l8).  The  calculation  of  energy  as  a  function  of  S  may  he  carried 
out  in  complete  analogy  with  the  procedure  in  Section  III  B  above .  Thus 


Eq.  (3*0  must  be  replaced  by: 


ItM*  + 


and  Eqs.  (39)-(4l)  by: 


Ji(/  +/  )<K*iW i 

•  l-fc 


-(['*  + f“)  Wf-*)  U'UV  (54) 

(55) 

e  --  ijjl'f)  .  (56) 

In  analogy  with  (39a)-(4la),  it  is  possible  to  derive  an  equivalent  set  of 
equations  where  the  Integrals  run  from  -a  to  a  : 

Mi)  =  t.h)  +  (54.) 


(55a) 
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e,  =  r, 


When  a  is  small,  we  have,  approximately, 


ft 

€. 


7  * 


/V/  JL.  — 


T  H* 


or 


€,  ^  ITp, 


(56a) 


(57) 


(58) 


provided  Is  small. 

The  asymptotic  values  of  t  and  p  when  a  is  large  may  be 
obtained  as  follows  using  (5*0 “(56).  First  rewrite  (5*0  in  the  form; 

f(V*  J/frJW*  (59) 

2 

Since  K(z)  decreases  as  l/z  for  large  values  of  z,  it  is  reasonable  to 
neglect  K(§  +  T|)  relative  to  K(§  -  T])  on  the  right  hand  side  of  (59)> 
provided  a  is  large  and  we  are  only  Interested  in  values  of  §  greater 
than  or.  Thus  for  large  or  we  have,  approximately, 

ft's)  *  UxV)'‘-fj((T-v)f(v)<ii. 

Let  §  s  ax,  T)  =  cry,  and  f(crx)  -  p(x).  Elation  (60)  becomes: 

r  t* 

)»(«)  J  *.K 


(61) 
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Now  as  a  becomes  very  large,  the  kernel  oK(ax  -  ay)  approaches 
a  Dirac  delta  function,  6(x  -  y).  If  ye  replace  It  with  a  delta  function 
inside  the  integral,  the  solution  of  (6l)  is: 


y(*)  *  (  71  d*  %*) 


-/ 


(62) 


Thus  p(x)  is  a  function  of  <y,  but  only  through  a  multiplicative  constant; 
otherwise  it  is  a  smoothly  varying  function  of  x.  This  provides  an  a 
posteriori  justification  for  replacing  the  kernel  in  (6l)  with  a  delta 
function.  Finally  we  obtain  for  p  and  e  the  results: 


jO  ~  i  /  (  x  *) 
«  ~  i/(3  r, 


and 


(63) 


(64) 


when 


is  small. 


^ 
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D.  Numerical  Solutions  of  the  Integral  Equations 

The  integral  equations  (39)  and  (54a)  were  approximated  by  a  set  of 
4l  coupled  linear  algebraic  equations,  and  these  vere  solved  on  a  high 
speed  digital  computer.  Ihe  results  are  shown  in  Fig.  1.  The  normalized 
energies  «  ,  measured  from  the  ferromagnetic  ground  state,  and  , 
measured  from  the  antiferromagnetic  ground  state,  are  plotted  as  a  function 
of  =  S/N,  where  S  is  the  total  spin  of  the  state,  and  p  *  \  -  • 

Note  that  =  .693  -  e-  The  lower  curve,  computed  by  means  of  (39)-(4l), 
represents  not  only  the  lowest  energy  for  a  given  p^  for  states  in 
class  C,  but  should  also,  according  to  the  discussion  In  Section  II,  be 
the  lowest  energy  for  sill  states  in  the  antiferromagnetic  chain.  The 
open  circles  In  Fig.  1  are  the  lowest  energies  for  an  antiferromagnetic 
chain  containing  9  spins,  while  the  solid  circles  are  the  corresponding 
energies  for  a  chain  containing  10  spins.  These  energies,  computed 
directly  from  the  Hamiltonian  (l)^  and  therefore  Independently  of  the 
Be the  formalism,  fall  surprisingly  close  to  the  solid  curve --which  should 
represent  the  limit  as  N  becomes  infinite --and  thus  provide  additional 
support  for  the  "theorem”  of  Section  II. 

The  numerical  calculations  provide  some  confirmation  for  the  asymptotic 
estimates  contained  in  (48),  (49),  (50),  and  (51)*  For  a  -  3*0,  the 
numerical  results  are: 

-  '  fsi  mo'*  »  .*0+  e'LCr< 

(65) 

6,  ~  4.5 7 pi  ,  .  (i 
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Equation  (52)  is  also  borne  out  by  numerical  calculations.  For  instance, 
at  a  =  .2: 


(>-  .  Ilf 2  -  .  92/  (2*/^) 


(66) 

*  =  *  /.57^ 

The  upper  curve  in  Fig.  1  represents  the  maximum  energy  for  states 
in  class  C  and  was  computed  by  means  of  (54a)-(;>6a).  The  calculations 
were  supplemented  by  the  use  of  (54)-(56)  for  larger  values  of  or.  It 
must  be  emphasized  that  this  curve  has  significance  only  for  the  class  C 
and  not  for  the  entire  energy  level  spectrum,  as  in  general  (except, 
perhaps,  for  p  very  close  to  0  or  £)  there  will  be  states  in  class  B 
lying  both  above  and  below  the  curve.  Once  again,  the  numerical  calcu¬ 
lations  confirm  the  asymptotic  behavior  for  small  and  large  a  values 
given  by  (57),  (58),  (63),  and  (64).  Thus  for  a  =  .2 
*s  .  OS’  I  4  S'  •»  t  .O  Z9  C*  /  4) 

e,  »  .!S9\  9H  (7;  fr)  (67) 

and  for  a  =  10. 


(?  -  .0  33  63  *  t.osy 

6-  Isn  x/o" 4-_ / .00 1  (2  nxp*/3). 


(68) 
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IV.  Itegnetic  Moment  and  Susceptibility  at  Zero  Temperature 

In  the  presence  of  a  magnetic  field  H  along  the  positive  z  axis, 
there  is,  in  addition  to  the  exchange  energy  (1),  a  Zeeman  energy : 

s  1mH  ?  Si*  (69) 

i 

where  g,  is  the  Bohr  magneton,  g  the  election  g  factor,  and  S^z  the  z 
component  of  spin  for  the  i'th  atom.  The  lowest  energy  level  of  the  chain 
for  a  given  value  of  total  qpln  S  in  the  presence  of  a  magnetic  field  will 
be: 


(S}=  £  (5)  -jmHS 

where  E(s)  is  the  lowest  energy  in  the  absence  of  a  magnetic  field, 
thus  equal  to: 

e(s)‘l«Te,ts/n) 

By  «1(S/H)  «  e^(p1 )  we  mean  the  function  corresponding  to  the  lower 
curve  in  Fig.  1  and  determined  by  Eqp.  (39a)-(4la).  Next,  let  us  choose 
among  sill  the  values  of  S  the  value  Sq  for  which  (70)  is  a  minimum.  This 
can  be  done  by  differentiating  the  right  hand  side  of  (70)  with  respect 
to  S  and  setting  the  derivative  equal  to  zero,  with  the  result: 

l  T  «/  li./n)  (72 


(TO) 

and 


(71) 


*  *» -vex** 


26 


At  zero  temperature  the  free  energy  of  the  chain,  F,  is  equal  to 
Ern1n(So).  From  this  one  can  calculate  the  average  magnetization  per 
spin : 

M-  df/dH  *  (73) 

The  magnetization  as  a  function  of  applied  field  is  plotted  in  Fig.  2. 

At  fields  in  excess  of  Uj/g^,  the  magnetization  is  completely  saturated. 

There  is  a  veil  defined  susceptibility  for  values  of  H  much  less  than  J/gu: 

K*  J1/** /(  T)  .  (7*0 

The  numerical  constant  appearing  on  the  right  Bide  of  (7*0  is  determined 
by  the  asymptotic  expression  (51),  and  is  only  known  to  a  precision  of 
about  one-half  percent,  as  indicated  in  the  discussion  preceding  (51) . 

Strictly  speaking,  the  result  (7*0  holds  only  for  the  infinite  anti¬ 
ferromagnetic  chain.  For  chains  containing  a  finite  number  of  atoms  the 
susceptibility  at  zero  temperature  in  the  limit  of  zero  magnetic  field  is 
either  zero  (for  N  even)  or  infinite  (for  H  odd).  Calculations  at  finite 
temperatures  for  short  chains  containing  up  to  11  spins'*"0  indicate  that 
(7*0  is  probably  not  an  unreasonable  limit  for  the  infinite  chain. 

Hulthrfn11  derived,  by  means  of  a  "statistical  model,"  an  approximate 

value  for  x  ,  identical  with  (7*0  except  that  the  constant  (2tt2)-1  *  .0506 

12 

is  replaced  by  .0593-  The  "self  consistent”  procedure  used  by  Thin  yields 

O  .1 

in  place  of  (2n  )  the  constant  .0556,  only  10£~ larger  than  the  exact  value. 
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Appendix  A.  Existence  of  Solutions  of  Egiation  (6)  for  States  In  (Hass  C 

We  shall  prove  that  the  transformation  (l4),  under  suitable  restrictions 
possesses  a  fixed  point,  and  hence  Eq.  (6)  a  solution.  Let  V  be  the  real 
r-dimensional  space  of  vectors  of  the  form  (k^,  k£,  . kr)  and  let  K  ' 
be  the  subset  of  V  consisting  of  those  vectors  vbose  components  satisfy 
the  inequalities: 

k,  £  ir,  V' 

£  tn  (l~  IV'*)  .  (Ai) 

-  kj  i  in  A/'» 

The  function  (see  (7))  is  positive  for  k^  >  ky  negative 

for  k  <  k  and  continuous  for  all  k.,k  in  the  interval  (0,2n)  except 
to  J  t 

at  the  points  k,  ■  k  .  Hence  the  transformation  (l4)  restricted  to  the 

V  t 

subset  K  is  continuous.  Furthermore  the  function  f(k.,k  )  in  either 

the  region  k  <  k  or  the  region  k.  >  k  is  monotone  increasing  in  the 
J  t  0  * 

first  argument  and  monotone  decreasing  in  the  second.  Utilizing  this  fact 
one  obtains  the  inequality: 

kh,  '  k/  *  *r  V  USh-  -in-1  .  m 

For  states  in  class  C  ve  have  the  ine quality  (13)  which  together  with  the 
fact  that  j<  tt  implies  that 

hilt  *  IHA/'I  (A3) 

The  aforementioned  properties  of  f  together  with  condition  (11)  guarantee 


fir*----  *» 


the  inequalities: 


V  1  In  V"' 

kr'  UsO-V'J 
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(A4) 


A  comparison  of  (A3)  and  (AU)  with  (Al)  shows  that  the  transfonaation 

(lU)  carries  the  set  K  into  itself.  The  set  K  is  closed,  convex,  and 

bounded;  this  together  with  the  continuity  of  (14)  on  the  set  K  implies 

1^ 

the  existence  of  at  least  one  fixed  point  by  Brouwer's  theorem. 
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Appendix  B.  The  Ledinegg-Urban  Calculation 

For  convenience,  we  shall  enclose  equation  numbers  from  the  paper 

k 

by  Ledinegg  and  Urban  in  square  brackets,  and  transcribe  the  equations 
in  our  notation .  Their  paper  is  concerned  with  the  problem  of  obtaining 
asymptotic  values  for  p  and  c  (or  p1  and  *1)  as  a  function  of  a  in 
the  limit  of  large  a,  using  (39)-(4l) •  Their  results  are: 

~  .  fIS  6  •*’  **  [38d] 

**  ^  / it1  [37] 

The  result  [38d]  is  not  much  different  from  our  calculations — see  (1*8) 
and  (50)— except  we  obtain  a  smaller  value  for  the  numerical  constant. 
However,  a  comparison  of  [37]  with  (49)  shows  a  very  important  difference: 
The  exponent  which  appears  in  the  latter  is  twice  the  value  in  the  former, 
and  the  factor  a~  is  missing.  The  numerical  solutions  to  Eqs.  (39)-(4l) 
support  (49)  rather  than  [37]. 

Ledinegg  and  Urban  obtain  their  asymptotic  estimate  as  follows.  The 
function  f(§)  in  (39)  (their  Eq.  [23])  is  set  equal  to: 

*(V  *  /.  I'S)  +■  W*)  [2*a] 

where  fQ(§)  is  the  solution  of  (24)  given  by  (32).  Hext  they  assume  that 

I WD I  /  4  (*)  «  1/  [2‘o;i 
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stating  that  this  inequality  will  be  justified  through  an  explicit  calcu¬ 
lation  of  the  function  It  is  not  stated  whether  the  inequality  [24c] 
is  to  hold  for  |g|<  a  or  for  |g|  >  or.  They  make  explicit  use  of  [24c] 
for  |gj>  a  in  deducing  the  approximate  integral  equation  for 


*(1)  *  (  [*  +/i<  )  4  (i-v)  t.l«) 


[26b] 


Their  justification  for  [24c],  however,  is  contained  in  the  equation: 


Ws)  ««*/'•  til*--*)]  [3»] 

where  $  is  a  relatively  slowly  varying  function  of  its  arguments. 

This  shows,  indeed,  that  the  inequality  [24c]  is  Justified— but  only 
for  |§|  <  cr  .  It  is  easy  to  prove,  using  their  expression  [33d]  for  $, 
that  for  |§|  »  a,  [24c]  is  not  only  incorrect,  but  the  reverse  inequality 
holds,  since  for  large  g  , 

K  (J)  ~  e  *  “  ^  (44) 


However,  it  is  for  the  region  |g|  >  cr  that  the  inequality  [24c]  must 
be  invoked  in  order  to  Justify  [26b].  Thus  we  infer  that  their  treatment 
is  not  Internally  consistent. 

It  is,  in  fact,  very  difficult  to  attack  (39)  by  a  straightforward 
perturbation  approach,  since  which  is  of  the  order  of  exp(-rrcr),  turns 
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out  to  be  the  difference  between  two  much  larger  terms,  both  of  order 
exp(-^rr<*).  The  advantage  of  (39a)-(4la)  over  (39)-(4l)  for  the  case  of 
large  or  consists  precisely  in  the  fact  that  no  such  near  cancellation 
occurs  In  calculating  either  *1  or  q^.  Even  the  zeroeth  order  approxi¬ 
mation  to  (39a) --replacing  f($)  by  fQ(£)--leads  to  the  correct  functional 
form  for  the  asymptotic  values  (48)  and  (49),  although  not  to  the  correct 
numerical  coefficients. 


Figure  Captions 


Fig.  1.  Highest  and  lowest  energies  for  states  In  class  C  as  a  function 
of  the  total  spin  S.  The  circles  represent  the  lowest  energies  for 
finite  chains  containing  9  And  10  atoms. 

Fig.  2.  The  magnetization  as  a  function  of  field  for  the  antlferro- 
roagnetlc  linear  chain  at  zero  temperature. 
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